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y-^ , We describe the interplay of quantum and thermal fluctuations in the inflnite-range Heisenberg 

("■^ ■ spin glass. This model is generalized to SU{N) symmetry, and we describe the phase diagram 

fvj ' as a function of the spin S and the temperature T. The model is solved in the large A*' limit 

and certain universal critical properties are shown to hold to all orders in 1/N. For large 
S, the ground state is a spin glass, but quantum effects are crucial in determining the low 
T thermodynamics: we find a speciflc heat linear in T and a local spectral density of spin 
excitations, x'l'od^) ~ ^ ^^r a spin glass state which is marginally stable to fluctuations in 
>D ' the replicon modes. For small S, the spin-glass order is fragile, and a spin-liquid state with 

x'loc ~ tanh(a;/2r) dominates the properties over a significant range of T and lj. We argue 
that the latter state may be relevant in understanding the properties of strongly-disordered 
transition metal and rare earth compounds. 

i . . 

C/3 ' The study of intermetallic compounds of the transition metals and rare earths has been a subject at the forefront 

^ ; of condensed matter physics for some time now [^ || . A rich and complex variety of behaviors is observed in low 
3 ' temperature electrical and magnetic measurements, much of which lacks a comprehensive theoretical description. The 
complexity arises from the dominant role played by the local magnetic moments on the d and / orbitals and their 
interactions with each other and the itinerant charge carriers. 

It is convenient to begin our discussion in a phase with well-established magnetic order, in which each magnetic 

Q ■ moment is effectively static. This static moment could be polarized in a regular manner (as in a commensurate 

CJ ' antiferromagnet or an incommensurate spin density wave), or point in random directions (as in a spin glass state). 

In most realistic systems, the magnetic moment is either quite small, or has averaged to zero by dynamic quantum 

fluctuations: so it is useful to consider mechanisms which reduce the magnetic moment, and eventually cause it to 

^ , vanish at a quantum phase transition to some paramagnetic state. Two distinct routes to such a quantum phase 

00 ■ transition can be envisaged, and, we believe, the interplay between them is at the heart of the complexity of the 

OO [ problem. In the first route, originally discussed by Doniach ||], the moment is quenched by Kondo screening by the 

CO itinerant electrons: theories of such quantum critical points have been proposed] |, |5|, |g, |7| in which the predominant 
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role of the itinerancy is to overdamp the collective magnetic excitations. In the second route, the exchange interactions 
between the moments play a more fundamental role: a pair of spins interacting with an antiferromagnetic exchange 
^-^ ' prefers to form a singlet valence bond, and the proliferation of such singlets can destroy the magnetic order. Analytic 
■"■^ — , theories for such transitions has been made mainly for systems without quenched disorder]^. Simple models of 

d • crossovers between these two routes have also been presented |^, |l^ |ll|. 

^ ' This paper will present a detailed study of the second route to destruction of magnetic order for the case of a 

strongly random system with spin glass magnetic order. There are a number of motivations for focusing on random 
' O ' systems. First, randomness is inevitably present in all materials, and it is clear that it strongly perturbs the low 

C , temperature properties. Spin glass order is present in a number of systems, while others appear to be in the vicinity 
of such a state. Finally, a technical motivation is in the structure of the mean-field theory we shall present: it builds 
in important feedback effect between the inter-site magnetic correlations and the single-site spin dynamics, and this 
is crucial to all the non-trivial spin correlations we shall describe. Such a feedback is absent in previous studies of 
the magnetic quantum critical point, and it has been argued that this is an important limitation for them[[7|, |ll|]. A 
different route to incorporating these feedback effects has been taken in some recent studies|l2|; however, they discuss 

C^ only the paramagnetic state of their model, and the extent to which magnetically ordered states preempt their results 
remains to be clarified. 

This paper is organized as follows : In Section |, we present our spin glass model and give an outline of our results, 
including the phase diagram. Section O is devoted to the nature of the paramagnetic solutions, and more specifically 
to the quantum critical regime. Section |l| is devoted to the spin glass phase and to the various regimes within this 
phase, as a function of temperature and of the size of the spin. The Appendices contain technical details and some 
additional results on the quantum rotor and Ising spin glasses of Ref. O. 



I. MODEL AND OUTLINE OF THE RESULTS 

The numerous recent studies of quantum fluctuations in spin glasses p4| , have focused either on infinite-range 
models of Ising and rotor models [O, Eq| or models in low dimensions which flow to strong disorder fixed points 
p6| , O, ll^. Here we shall continue the study of infinite range models, but will consider a model of Heisenberg spins 
: in this case, the path integral for each spin has a important Berry phase term which imposes the spin commutation 
relations. As we will see, this leads to a great deal of new physics|l9| and non-trivial dynamic spin correlations even 
in the spin glass state. More specifically, we present a complete solution of the quantum Heisenberg spin glass on 
a fully connected lattice of Af sites with strong Gaussian disorder, both in the paramagnetic and the glassy phase, 
when the spin symmetry group is extended from SU{2) to SU{N) and the large-A^ limit is taken. In the limit of 
large connectivity, (dynamical) mean-field techniques apply and the model can be reduced to the study of a self- 
consistent single-site problem, which is however still highly non-trivial because of quantum effects. The large- A^ limit 
is instrumental in allowing for an explicit solution. Nevertheless some of our results regarding the quantum critical 
regime have been extended beyond the large-N limit. In a recent publication [EOJ, we summarized the main results of 
the present study. Here, we provide detailed derivations and new results, such as a full discussion of the paramagnetic 
phases and a discussion beyond large- iV. 

The model considered in this paper is defined by the Hamiltonian: 

^-^^T.J^^^^-^j^ (1) 

where the magnetic exchange couplings Jy are independent, quenched random variables distributed according to a 
Gaussian distribution 



As already pointed out by Bray and Moore [^, after using the replica trick to average over the disorder [g2|, the 
mean-field (infinite dimensional) limit maps the model onto a self-consistent single site model with the action (in 
imaginary time r, with /3 the inverse temperature) : 

Seff = Sb-^J^ drdr' Q'^\t - r') 5^^^) ' ^\r') (3) 

and the self-consistency condition 
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where a,b = 1, • ■ • , ri denote the replica indices (the limit n — > has to be taken later) and Sb is the Berry phase 
of the spin [|l9|. Due to their time-dependence, the solution of these mean-field equations remains a very difficult 
problem for N — 2, even in the paramagnetic phase. Thus, in Ref. pl| as well as in most subsequent work p3[ , 
the static approximation was used, neglecting the r— dependence of Q'^t). This approximation may be reasonable 
in some regimes but prevents a study of the quantum equilibrium dynamics, and is particularly inappropriate in 
the quantum-critical regime. However, this imaginary time dynamics has been explicitly studied in a Quantum 
Monte Carlo simulation in the paramagnetic phase with spin S = 1/2 by Grempel and Rozenberg pJ. Recently, 
we introduced a large- A^ solution of the mean field problem ||20|] , in which the problem is exactly solvable and, as 
explained below, the solution provides a good description of the physics of the N — 2 mean field model, to the extent 
the latter is understood. More specifically, in the following we will consider two different types of spin representations 
for the SU{N) spins : 

a) Bosonic representations : the spin operator S is represented using Schwinger bosons b by Sap = bl^bis — SSaf3, 
with the constraint X^a ^q^" — SN (0 < 5"). In the language of Young tableaux, these representations are 
described by one line of length SN. They are a natural generalisation of an SU{2) spin of size S. 

b) Fermionic representations : the spin operator S is represented using Abrikosov fermions / by Sap — fafp^Qo^ap, 
with the constraint J^a fafa = QqN {0 < qo < 1). In the language of Young tableaux, these representations are 
described by one column of length qoN. Note that for SU{2), only S = and S = 1/2 can be represented in 
this manner. 



In the following, we refer to the model with bosonic (resp. fermionic) representations as the bosonic (resp. fermionic) 
model. In the fermionic model, quantum fluctuations are so strong in large- iV that the spin glass ordering is destroyed 
[ fi9[ , contrary to the bosonic model, where a spin glass phase exists, as explained below. The two models have 
different theoretical interest : if one wants to concentrate on the quantum critical regime, above the spin glass 
ordering temperature, one can use the fermionic model (as e.g. in Ref. Eq). However, since we are interested in the 
spin glass phase itself, we will now focus on the bosonic model. Nevertheless, our results on the paramagnet will be 
valid for both cases with only slight modifications explicitly quoted below. 

In the N —> oo limit, the mean field self-consistent model (|3|) reduces to an integral equation for the Green's function 
of the boson Gf'{T) = —{Th°-{T)W^{Q)) where the bar denotes the average over disorder and the brackets the thermal 
average [|9J : 

{G-^r\iUr,) = iv^Sab + X'^Sab - Sf (i;/„) (5a) 

Ef (r) = j'{Gfir)fGf{-T) (5b) 

Gr{T = 0-) = -S (5c) 

Similarly for the fermionic model, we have : 

iGj'r\lU;n) = iiOnSab + X^Sab ~ T.f {lUn) (6a) 

Ef (t) = -J^{Gf{T)fGf{-T) (6b) 

Gf(T = 0-)=go (6c) 

In these equations, v (resp. lo) are the bosonic (resp. fermionic) Matsubara frequencies and the inversion should be 
taken with respect to the replica indices a, b. Note that our conventions for the sign of the Green functions in this 
paper slightly differ from those of Ref. O. Note that, although the equations are written in term of G, the physical 
quantity is the local spin susceptibility xiociT) = {S{t)S{0)) which is given in the large- A^ limit by 

Xiocir) = GriT)Gri-r) (7) 

In many instances below (where we mostly focus on the bosonic case), we shall drop the index 6 in Gb- 

From both analytical and numerical analysis of these integral equations, we have constructed the phase diagram 
displayed on Figure 0, as a function of the size of the spin S and the temperature T. Let us give here a brief overview 
of the main features of this phase diagram, which will be studied in great detail in the rest of this paper. As is evident 
from Fig H^, it is useful to divide the discussion into models with S large and S small. Both regimes are accessible 
in the large N limit, where S is effectively a continuous parameter taking all positive values. For the physical case 
A^ = 2, we will present evidence later that at least S — 1/2 is in the small S regime for the infinite-range model; 
moreover, we can expect that at least some of the consequences of increased quantum fluctuations in a realistic model 
with finite-range interactions are mimicked by taking small S values in the large N theory of the infinite-range model. 
For large S, the ground state must clearly be a spin glass (Fig[l]). However, even for very large S, it is necessary 
to consider quantum effects in understanding the lo w T excitations and thermodynamics, and these have not been 
previously described. In this paper (Sections |IIIC2 and [II D| ) we will show that the local spin susceptibility has a 



low energy density of states which increases linearly with energy. At the same time, the specific heat also has a linear 
dependence upon temperature. These results hold for temperatures T < J\/S, although characteristic excitations 
have an energy of order JS for T < JS; we will provide scaling functions which determine the dynamic response 
functions at these energies. At even higher T, there is a phase transition to a paramagnet at T ~ JS^. For large S 
the static properties of this phase transition are well-described by a purely classical theory in which the 5^ in (m are 
commuting vectors of length S. Notice also that we indicate two critical temperatures, T^g and T^^: as we discuss in 



Section IH , these are a consequence of peculiarities in the nature of replica symmetry breaking, where the dynamic 
freezing into the spin glass phase {T^„), happens at a slightly higher temperature than the equilibrium transition 

For small S, we also find a spin glass phase (Fig P at T = 0, but the order vanishes at a small T. Moreover, its 
excitations and finite T properties are very different from those at large S. These are now dominated by signals of 



a "spin-liquid" state discussed in Ref. [19| and described in Section |II A| . In particular, we describe a novel quantum- 
critical region in the paramagnet where max(a-', T) is the characteristic energy scale, and the local dynamic spin 
susceptibility obeys Xiod^) ^ tanh(w/2T). We believe that aspects of this regime may be relevant to disordered 
transition metal and rare earth compounds in regimes where exchange interactions between the magnetic moments 
are playing a dominant role. Completion of this picture requires an understanding of the stability of the "spin-liquid" 
picture to mobile charge carriers, and this has also been addressed in a previous work ]25|. 
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FIG. 1: Phase diagram of the mean held bosonic model. There is a spin glass phase below the spin glass temperature 



T^ , which is determined with the marginality condition (See Section HI). T^^j is the spin glass temperature as 
determined with the stationarity criterion. 



II. PARAMAGNETIC PHASE 

Contrary to the classical case, the paramagnetic phase of quantum spin-glass models is non trivial in mean field 
theory. An early discussion of these solutions has been given in Ref. HO^, but we present here a much more complete 
description, and compare our results to the N = 2 case, when numerical results are available. Since in this section 
we look for paramagnetic solutions, we will consider only replica diagonal solutions of (0) : G"'' oc Sab- Two types 
of paramagnetic solutions have been found, that we will now consider successively : the spin-liquid solutions and the 
local moment solutions. 



A Spin-liquid solutions 

1 The large-N limit 

A low- frequency, long-time analysis of the integral equations (ph reveals that, under the condition that A — S(iO''') 
vanishes at low temperature, a solution can be found which displays a power law decay of the Green's function at 
long time pS( |: G{t) ^ Xj ^/t. These solutions display a singularity in the complex plane of frequencies, z, at z = 0, 
with an amplitude which can be parameterized by an angle as: 



G{z) 



j^^-iTr/4-i0 



for z ^ 0, Imz > 



(8) 



(Values of z on the imaginary frequency axis at the Matsubara frequencies will be denoted by i^„, while on the real axis 
will be denoted, uj.) Thus, these solutions display a slow local spin dynamics : Imxioci^) oc sgnuj for lu —^ 0,T — 0. 



Moreover, the local susceptibility xioc{T) = /p XIoc{t) dr diverges as xiociT) ~ InT/J at low temperature. More 



precisely, one can find the thermal scaling function characterizing the r 
paper |E5| : 



oo, T — > limit, as explained in a previous 



XIoc{t,(3) cx 



n/P 



sin TIT / (3 



Jx'i'oci^,T) cxtanh — 



(9) 



Note that in the paramagnetic phase of quantum models, the local susceptibility Xioc{T) (which is the response to a 
local magnetic field) differs from the uniform susceptibility x{T) (response to a constant magnetic field), contrary to 
classical spin glass models where x = Xioc p^ : this is a consequence of the commutation relations of the spin, as can 
be seen for example in the high-temperature expansion in the SU{2) model. In this large- A^ limit, it can be shown 
that xiT) ^ Xioc{T) for T ^ and numerical computations indeed suggest that x{T) ^ const. p5[ . 

Remarkably, the parameter 9 which characterizes the spectral asymmetry |2(| of the spectral density at low frequency 
can be explicitly related to the size S of the spin (which involves a priori an integral of the spectral density over all 
frequencies). This is very similar to a kind of Friedel sum rule applying to this problem, and indeed the derivation 
follows a very similar route, based on the existence of a Luttinger-Ward functional. (Interestingly enough, the 
"boundary term" which usually vanishes in such derivations contributes here a finite value). This derivation is 
presented in detail in Appendix Rl where the following relation between 9 and S is established: 



9 sin 29 

n 4 
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in the bosonic model 
in the fcrmionic model 



(10) 



This relation has important consequences for the physical properties of the spin-liquid solutions. First, we note that 
the spectral density must obey the positivity conditions: ImG'/(w -I- iO"*") < and sgn{uj)ImGb{uJ + iO"*") < 0. Hence, 
in the fcrmionic case, 9 must obey —j<9<j. It is easily checked from (noh that 9 precisely describes this range 
of parameters as qa is varied from qo — to qo = 1, and that the 0(qo) relation is unique. This suggests that the 
spin-liquid solution is an acceptable low-temperature solutions for the whole range of qo in the fcrmionic case. In 
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FIG. 2: S as a function of 9 for the bosonic model. The solid line is given by relation ( |JC|) and the points were 
obtained previously from a numerical solution of the saddle-point equation at zero temperature nj 

contrast, in the bosonic case, the plot in Fig.|| shows that ( [lOI) actually defines two values of 9 (in the allowed range 
J < < ^) for a given spin S as long as 5 < Smax — 0.052, while no value of 9 is found for S > Smax- This implies 
that no paramagnetic solution of the spin-liquid type is found at zero-temperature in the bosonic case as soon as 
S > Smax (note that furthermore Smax is very small). For S < Smax, such solutions exist at zero temperature (even 
though they are not the true ground-state, see below) with the locally stable solution corresponding to the smallest 
of the two values of 9. However, even for S > Smax, at low (T < J) but finite temperature, the spin liquid solutions 
do exist in the bosonic model. By this, we mean that a numerical computation in imaginary time gives a solution 
which exhibits the scaling form (O), for which a unambiguous value of the spectral asymmetry 9 can be defined and 
computed numerically. At very low temperature, these solutions are unstable to the spin glass solution, but above 



the spin glass temperature at low spin, they are relevant in the quantum critical regime associated with the quantum 
critical point at 5" = 0. We shall comment in more detail, at the end of the following section, on the nature of the 
paramagnetic solutions found at low temperature for small values of S, in the bosonic case. 

Another consequence of relation (10) is that it allows to predict that these spin- liquid solutions have a non-zero 
extensive entropy at zero-temperature and to calculate the value of this entropy analytically. The derivation of this 
result follows very closely a similar analysis of the overscreened multichannel Kondo problem in the large-N limit, 
performed in Ref. Eq, B^ and only the main steps will be repeated here. This can be done either in the bosonic model 
or in the fermionic one, with slight modifications. Since the spin liquid solutions are relevant at zero temperature 
only in the fermionic model, we shall present the result in this case. First, denoting by S the value of the entropy per 
spin at zero temperature, one establishes the following thermodynamic equality: 



dS 
dqo 



dX. 



\T=0 



(11) 



Then, a low temperature expansion is used which allows to relate the slope of A(T) to the spectral asymmetry 
parameter above, so that one finally gets (in the fermionic case): 



dS _ sin(7r/4 - 6) 
dqo ^ ^ sin(6' + 7r/4) 



(12) 



The entropy is then obtained by integration over the size of the spin, with the physically obvious boundary conditions 
'5(<Zo = 0) = S{qQ = 1) = 0. The resulting value of the entropy as a function of qq is plotted in Fig.|3. 




FIG. 3: Entropy as a function of the size of the spin (qq) in the fermionic model. 



Finally, we comment on the physical nature of the spin-liquid paramagnetic solutions found in this section. These 
solutions correspond to a partial screening of the local moment at each site, due to the interaction with the other 
spins. As a result the local susceptibility diverges logarithmically (much slower than a Curie law), but an extensive 
entropy is still present at T = 0, indicating a degenerate state. From a local point of view, the physics is somewhat 
similar to an overscreened Kondo system, but here the gapless bath which quenches the spin is not external but self- 
consistently generated by the other spins. We suspect that the physics of this phase has to do with the degeneracy 
of the (large- iV generalization) of th e "trip let" state in which two spins are bound whenever a strong ferromagnetic 



bond Jij is encountered. In Section II A 2, we show that this spin- liquid regime is not a peculiarity of the large- iV 



limit but indeed survives in the mean-field description of the quantum critical regime of a SU{2) quantum Heisenberg 
spin-glass. It would be very valuable to gain a more direct understanding of this gapless spin-liquid regime from 
a study of the problem for a fixed configuration of bonds, before averaging over disorder. This could be achieved 
numerically and is left for future studies. 



2. Beyond the large-N limit 

This subsection will show how recent renormalization group analysis of related models |0, Eq, 29 imply that the 
above spin-liquid solution applies to all orders in 1/N . In particular, the large A'^ solution with Ivaxioc oc sgn uj 
for small oj acquires no corrections to its functional form: the only changes are to the non-universal proportionality 
constant. All the discussion below will be in a paramagnetic phase where it is sufficient to consider only a single 
replica, and so we will drop replica indices in this subsection. 



We begin by rewriting (p) in the following formjll 



Seff = Sb-1o f dT^ir) ■ ^(r) (13) 



where 70 is a coupling constant and cj) is an annealed Gaussian random field with ( iji (t) ■ cj) (0)) = 1/|tP '^. It is 
reasonable to expect that the spin correlations in the quantum ensemble defined by ( |l3|) decay with the power-law 
{S{t) ■ 0(0)) ^ l/|r|'^, and we are interested in determining the value of the exponent a. A simple extensionlQ of 
the solution discussed above implies that in the large N limit a = e. Here we will argue that this equality is in fact 
exact for all N. Now using the self-consistency condition ( |4|) w e obtain a = e = 1, which then implies Imx;oc oc sgn uj. 
The field-theoretic renormalization group analysis of ( |13| ) was discussed in Ref. ES, and we will highlight the 
main results. The key observation is that renormalization of the theory ( |l3|) requires only a single wave-function 
renormalization factor Z , and that there is no independent renormalization of the coupling constant 70. This result 
was established diagrammatically in Ref. 09, and we will not reproduce the argument here. So if we renormalize 
the spin hy S — vZSr, then the coupling constant renormalization is simply 70 ~ /i''^7/vZ, where /x is a 
renormalization scale. The renormalization constant is in general a complicated function of 7, and was determined to 
two-loop order in Ref. |2^: 

272 74 
Z=l '- + ^ + ... (14) 

in a minimal subtraction scheme. However, even though Z is not known exactly, the exponent a can be determined 
exactly. Standard field-theoretical technology shows that the above renormalizations imply the /3-function 

«.)^-?ri4i^v'^ (15) 



2 V 2 91n7 
Furthermore, the exponent a is given by the value of 

<t(7)-/3(7)^ (16) 

at the fixed point 7 = 7* where (3{'j) vanishes. Comparing dl3) and (fiq) we see that 

/3(7) = -(e-a(7))7/2. (17) 

Clearly, a zero of the f3 function must have a = e, and this establishes the required result. 

We note that similar examples of a critical exponent being valid to all orders (in spite of a non-trivial /3-function) 
can be found for other models in the statistical mechanics of disordered systems (see e.g |^, ^ ) . 

B "Local moment" solutions 

In a mean field model, one usually expects to find locally stable (while possibly unstable to ordering) paramagnetic 
solutions of the mean-field equations down to zero temperature. Hence, the absence of solutions of the spin-liquid 
type for S > Smax suggests that a different kind of paramagnetic solution should exist for those values of the spin. 
Indeed, we have found that the integral equations (|5|) have another class of paramagnetic solutions in the bosonic 
case. These solutions actually exist for all values of the spin S and down to zero temperature. Hence, they coexist 
at low S with the spin-liquid solutions in some range of temperature. Their physical nature is very different from 
the previous spin-liquid solutions, and as discussed below they are not very physical solutions when considered at 



8 
low temperature. They are characterized by a Green's function which does not decay at long times and obeys the 



asymptotic behavior: Gb{T) ~—S^ 



- e '^ '^•^ ■^. In contrast to the spin- liquid case, A diverges for T 
A 



T 



in this regime: 

(18) 



Finding numerically these solutions of dq) requires some care. We have used an algorithm in which we solve (@) 
in imaginary time for G{t), for a fixed value of r = J„ G{T)dT, and then adjust the number of particles to S' by a 
dichotomy on r. The local susceptibility xioc{t) obtained in this manner is displayed on Figure M. At high temperature, 




FIG. 4: XIoc{t) extracted from a numerical solution of saddle point equation (^ in imaginary time, for S* = 1, J = 1. 
The solid curve is low temperature (Jf3 = lOj, the dashed curve is high temperature (J(3 — .1); for intermediate 
temperatures, the curves interpolate between the two. 

we find xioc{T) — -^ — - as expected since the spin is essentially free. At low temperature, we find another Curie 
law, with a reduction of the Curie constant due to quantum fluctuations : 



Xioc(T) 



T 



for T -> 



(19) 



Hence, for these solutions, the effect of the interactions with the other spins is not strong enough to result in a 
qualitatively different screening regime, resulting merely in a reduction of the Curie constant. This is analogous to 
an under screened Kondo regime. 

These solutions are of a similar type than those found in in a quantum Monte-Carlo simulation of the SU{2) model 
Q. There also, a reduction of the Curie constant from S{S + l)/3 to 5^/3 was clearly observed. Thus, contrary 
to one of the conclusions of Ref. ^, the large- iV limit correctly reproduces the paramagnetic local moment solution 
found for the physical N ^ 2 case. Moreover, a numerical solution of the large-N integral equations ^ for real 
frequencies can also be obtained both at high and low temperature. The results are presented in Figure^ for both 
Pb(u]) = — l/7rImG'f,(u;) and Xiod'^) — ^ / ''^^''^Xioci^) ■ At high temperature, pb is centered around A r^ Tln((S'-f l)/5') 
and x'l'oci'^) i^ ^ simple peak. At low temperature, we find in x'l'od'^)-^ delta- function peak at zero frequency with 
weight S"^, which is associated with longitudinal relaxation (see Ref. gj for a discussion for N = 2) and two peaks, 
centered around ±A oc 1/T (Cf Eq. (Iq)) with a constant width, which are associated with transverse relaxation. 
Again, these results are very similar to the conclusions reached in Ref. ^ from a fit of the imaginary-time data for 
N = 2, using a model x {^)- The only difference is that the central peak is not broadened by thermal fluctuations in 
the large- A^ limit. 

Finally, since these solutions describe the formation of a local moment at low temperature, and that the onset of 
the quenching temperature (at which the reduction of the Curie constant sets in) turns out to be lower than the 
temperature where spin-glass ordering occurs (as shown in the next section), we consider these solutions as a mean- 
field artefact of the spin glass ordering. Static limits of such solutions are actually known to occur in the classical SK 
model. We also note that the internal energy of these solutions have an unphysical divergence as T ^ 0. 
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FIG. 5: Spectral densities for Gb (left) and for the local susceptibility xioc (I'iglit), for high temperature (top) and a 
very low temperature (bottom). These results are extracted from a numerical solution of saddle point equation ([gj 
in real frequencies. 

We close this section by noting that, for small values of the spin S, a rather intricate pattern emerges for the stability 
and coexistence of the two kinds of paramagnetic solutions described above. We have studied this numerically in some 
detail but do not report this here, since most of these phenomena occur below the spin-glass ordering temperature 
anyway. The important features have been displayed on Fig.n^. We emphasize again that the spin-liquid solutions 
are the relevant solutions describing the whole quantum-critical regime. Even though they are unstable at T = for 
S > Smax they remain consistent finite temperature solutions in a much wider range of values of S for T < J. At 
higher S and above the spin-glass ordering temperature, the paramagnetic phase behaves as a local moment, with a 
Curie constant getting gradually reduced as T is lowered. 



III. THE SPIN-GLASS PHASE 



In this section, we investigate spin-glass ordering in this model. The first observation that we make (Sec. Ill A) is 
that the spin-glass susceptibility (i.e the response to a spin-glass ordering field) is actually of order 1/A^ in the large-N 
limit. This does not preclude a spin-glass phase, but means that the transition is not associated associated with a 
linear instability. Indeed, we shall find explicit solutions in the ordered phase in the bosonic case, while the fermionic 
case does not have a spin-glass phase at iV = cx) (but does order as soon as 1/A^ corrections are considered). 



A Spin-glass susceptibility 



Here, we derive an exact expression for the spin-glass susceptibility, valid for arbitrary N in this mean-field model. 
The derivation is most conveniently performed using replicas. We note that in the presence of spin-glass order, the 
correlation function ( d "(t) • o ^(r') > acquires a non-zero value for a ^ b, which is however static, i.e independent 
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of r — r' (see e.g Ref. g). This crucial point is due to the fact that different rephcas are independent of one another 
before averaging, so that for a^h: 

{-^-{t) ■ ^"{r')) = (^^{r)) ■ {-^^{r')) = ('^'^{0)) ■ (^^0)) (20) 

In the foUowing, we shaU denote by qat the (normahzed) off-diagonal correlation function which is an order parameter 
for the spin-glass phase: 

9''^^]^E(^"W-^''(^')) (21) 

We consider the stability of the paramagnetic phase to this type of ordering, and introduce an ordering field Hat 
conjugate to qab- This has two effects: 

i) It adds to the effective action (ph an explicit term: 

5S=^ j dT j dT'Y,HabSa{T)-Sb{T') (22) 

a^b 

The normalization of Hab has been chosen in such a way that the change in the total free-energy is of order N . 

ii) It modifies the value of the self-consistent field Q""^ . The change of the off-diagonal component, to linear order, 
is imposed by the self-consistency condition (0) to be: 6Q"^'^ = Sqab, with 6qab the induced order parameter. 

One can then perform an expansion of the off-diagonal correlation function up to linear order in both Hab and 
SQ^K This yields: 

hab = ^xL {Hab + J'SQ'^') (23) 

in which xioc is the local susceptibility of the paramagnetic phase. Since tJQ"'' = Sqab, this finally yields the suscepti- 
bility to spin-glass ordering: 

^^_1 xL .24) 

^''-Hab Nl~{JxiocY/N ^ ' 

This formula has two important consequences. The susceptibility to spin-glass ordering is of order 1/A^ and any 
spin-glass instability at A'^ = oo must be associated with a non-linear effect of higher order (i.e come from terms of 
higher than quadratic order in qab in the free energy). Furthermore, (04) shows that for finite N, a (linear) instability 
into a spin-glass phase will occur when Jxioc{T) — VN. Hence the fermionic model, for which xioc diverges at low-T, 
will have a spin-glass instability for arbitrary large but finite N. More precisely, since the low-temperature behavior 
Jxioc ~ In ^ has been shown above to hold for all N for the paramagnetic solution of the fermionic case, we conclude 

from (|2J) that the spin-glass transition temperature depends on N in that case as : T/ ~ Je~^^ . 

B Spin glass solutions and the "replicon" problem 

We now turn to the explicit construction of solutions of the integral equations (^ with spin-glass ordering, in the 
bosonic case. The same reasoning as above shows that the Green's function G°'^{t) does not depend on r for a^h. 
Thus the most general Ansatz for the Green's function G°^ can be written as : 

C..(^)^fG(r)-5i (a = fe) (25) 

\ -gab (a 7^ b) 

where G{t) is a function of the imaginary time, gab a constant matrix and gi a constant. By definition, gi is fixed 
so that G is regular at T = 0, i.e. G{t) — > as t ^ oo. In the following discussion, we will restrict ourselves to 
solutions given by a Parisi Ansatz for gab (a replica symmetry breaking scheme). In the n — > limit (where n is the 
number of replica), this matrix becomes a function g{u) of a continuous variable u with < u < 1 123]. 
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Our equations (g) involve the Green function G, whereas the physical quantity is the susceptibility x"''(t) = 
G°'''{t)G°'''{~t). The order parameter qat widely introduced in the spin glass literature |g2[ is given here by Qab = <?a(,, 
or in the limit n — > : 

q{u) = g{uf (26) 

The Edwards- Anderson parameter is qEA — 9(1) — giXY [|l ^- Since at zero temperature, in the long time limit, 
we find lim^^oo G""(t)G""'(— t) — qEA we find gi = g{l), by definition of gi. More precisely, we look for solutions 
in which these two definitions of qEA coincide but it is not really an assumption in our computation : if this relation 
was violated, we would simply find for G a non vanishing limit for t — > cxd. 

Among the various possible replica symmetry breaking schemes ^^, we will now focus on one-step solutions, since 
we have not found any other, either two-steps or with continuous replica symmetry breaking. In this case, the function 
g{u) is piecewise constant : g{u) = g for < u < x and g{u) — g for x < u < 1. In the following, we will refer to x as 
the breakpoint. According to (^, the self-energy S has the form : 

I -J gL (a ^ b) 

with S is given by ( |31b| ) . Using the standard formulas to invert the Parisi matrices in the limit n — > p2| , we find : 

G-\lUn)^l>^n+X-i:{liyr.) (28) 

J^g^Giii^n = 0)^ = 1 + J^(3xg^G{iiyn = 0) (29) 

At this stage, it is useful to introduce here a new parameter O defined by : 

G{tJ^^)=--^ (30) 

Jg 



We then eliminate A in ( P8[ ) and G(0) in (29) and obtain finally a closed set of equations for G and g : 



[Gitl^n])' 



:^-(S(^^„)-I](0)) (31a) 



I](t) = J2 (g^{t)G{-t) - 2gG{T)G{-T) - gG^ir) + 2g^G{T) + g^Gi-r)) (31b) 

C{r = 0-) = -{S~g) (31c) 



-(- 



/3a^ = ^T ( S - Q ) (31d) 



The crucial observation at this point is that these saddle-point equations possess a one parameter family of solutions, 
parametrised by or equivalently by the break poin t x. This phenomenon already occur in other models which have 



a one-step replica symmetry breaking solution |33 . The determination of the breakpoint turns out to be the most 



difficult question of this analysis. Two possible criteria are : 

1. To minimize the free energy J-'{x) as a function of x, as would be required by the thermodynamics. In the 
following, we will refer to this as the equilibrium criterion. This criterion has been used in a previous attempt 
to understand this spin glass phase p3|. 

2. To impose a vanishing lowest eigenvalue of the fluctuation matrix in the replica space. We will refer to this as 
the marginality or replicon criterion. Although it is not really justified up to now, we will argue that it is the 
correct choice. 

This problem is not due to the quantum aspect of our model : it already appears similarly in some classical spin glass 
models, in the p-spin model for example. In this classical model, the study of the dynamics shows the existence of a 
dynamical transition T''^" above the static spin glass temperature T'"^ given by the static solution of the mean field 
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model. It turns out that, in this classical model, the replicon criterion has been proven to give the same transition 
temperature T"^ = T'^v"-, Moreover, it has been shown Q, ^ that the same phenomenon occurs in some quantum 
version of the p-spin model. Thus, using this condition, it is possible in some sense to mimic the dynamics by simply 
solving a static problem, although this is not fully understood at present. 

In the present model, the two criteria give a coherent solution but with totally different spectra of equilibrium 
fluctuations : the equilibrium criterion leads to a gap in x" (<£) whereas the replicon criterion is the only one which 
give a gapless x"{^) (^ similar observation was made in Ref. p3 in a one dimensional quantum model with disorder). 
We believe that in this quantum Heisenberg spin glass the replicon criterion provides us with the correct physical 
solution (Tc), contrary to the equilibrium solution, which gives the static transition temperature (Teg). However this 
claim cannot be proved in the present context : in particular, the static solution does give a full solution of (||). A 
study of the true Hamiltonian dynamics in real time and finite temperature of this quantum problem is necessary 
for a deeper understanding of this question, but this is beyond the scope of this paper. Let us now examine the two 
criteria separately in more details. 

1 The replicon criterion 

To apply the replicon criterion, we need to study the fluctuations of the free energy in the replica space around the 
one-step solution. In the large- A^ limit, the free energy is given by the expression : 

T[G-\ A] = i V Tr ln(zz.„ + A - ^-\ii^n)) +^y, f dr [G-\t)G'^\-t)\ ' - XS (32) 

Under infinitesimal variation 5gab for a ^ b^ the variation of the free energy is (up to second order) 

^^ ^'^ Mab^cdSgab^gcd- (33) 

a>b 
Od 

Strictly speaking, as this is a quantum problem, we have to simultaneously consider the variation of the diagonal 
component, SG{t) in ( p5| ) in a study of the fluctuations. In a spin glass phase, there is indeed a coupling between 
Sgab and SG{t) which modifies the fluctuation eigenvalues. Fortunately however, as we show in Appendix H, this 
coupling does not modify the eigenvalue ei and our main result (pq) below, and so we will neglect SG{t) here. The 
diagonalization of the n{n — l)/2 x n{n — l)/2 matrix M is briefly explained in Appendix H and gives three eigenvalues 

ei=3/3JV(l-3e2) 

e2 = ^^{e^~3 + 3(3Jg^e{l + e)) (34) 

63 = 6/3J252 [ipjg^Q - 1) 



A first consequence of this analysis is that replica symmetric solutions are unstable, since from Eq. (31d) they 
correspond to 6* = 1 and then ei < 0. Hence, these solutions will not be considered in the following discussion. 
A full solution of Eqs.(^) is required to show the positivity of 62, 63, but we immediately see that ei = for 

e« ^ -L (35) 



Quite remarkably, we will see below in Section [II C 2 and Appendix B 2 that precisely the same value of Q is 



selected by a criterion which is seemingly entirely independent. We will study the dynamic spectral functions in the 
spin-glass phase, as defined by G{t), and show that their associated spectral densities are non-zero as |u;| -^ only for 
the value of Q in (p3) . So marginal stability in replica space appears to be connected to a gapless quantum excitation 
spectrum. We may intuitively understand this as due to the availability of many low energy states when the system 
first freezes, but a better understand should emerge from a real-time analysis. 

2 The equilibrium criterion 

To apply the equilibrium criterion, we start from the expression of the free energy J-" and solve for x : 

^ = (36) 

ax 
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The computation of the total derivative ( |3q ) reduces to dx^{x)\^ ^ because the saddle-point equations (|^) for Unite 
n are equivalent to 



dJ" dF 



aC"'' dX 



= 



(37) 



as can be checked by an explicit calculation. Computing the logarithm in (|3^) (using Appendix II of Ref. |3^) and 
taking the derivative leads to : 



r'V-^K-.W... = o)) = --j^^ 



and finally to a equation for Q 



--i.4-?- 



(38) 



(39) 



This equation has two solutions : the replica symmetric one = 1 (unstable, as explained above), and a non trivial 
one O = Qeq ~ 0.4421 .... Contrary to the previous solution, we will see in Section IIIC2 that ImG has a gap for 
this value of 9. 



C The phase diagram 



Once Q has been determined, the equations ( pU| ) can be solved either numerically (both in imaginary time and in 
real frequency) or analytically in the S* — > oo limit. In the following, we will mainly restrict ourselves to = Qr since 
we believe that it is the correct solution. However all calculations have been redone for Q = Qeq with related results. 



1 Numerical solution 



First, the critical temperature, is obtained from the numerical solution of Eqs.(^^) in imaginary time : the spin glass 
order parameter q{T) = g'^{T) and the breakpoint x{T) are displayed in Figurelq as a function of the temperature. 
X increases linearly with T from at T = (there is no replica symmetry breaking at zero temperature) and Tsg is 




Tsg 



FIG. 6: The Edwards-Anderson parameter qea and the breakpoint x as a function of the temperature T for a fixed 
value of the size of the spin S* = 1 (J = 1). The transition to the paramagnet is given by the condition x = 1. 



determined by the condition x{Tsg) — 1, since we must have < x < 1 by definition |22]. Hence there is a discontinuity 
in q at the transition, but we will show below that the transition is second order. A careful numerical study shows 
that the transition is always driven by a; = 1 for all values of S and produces the critical temperature displayed 
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on Figure |[ The computation is similar for the critical temperature T^ given by the "replicon" criterion and for 
the critical temperature T^l^ given by the equilibrium criterion. Moreover, we find that T^„, the dynamic transition 
temperature, is higher than T^g, the static transition temperature : this is required by our physical interpretation 
of the two solutions but it was not obvious a priori from the integral equations solved. 

2 The large-S limit and spectral densities 

Further analytical insight into the spin glass phase itself can be obtained by considering various large-S" limits, 
which differ by the manner in which temperatures and frequencies are scaled with S (See Fig^). 

In a first simple large-S* limit, we take T large enough so that T / JS^ is of order unity. This is the simple classical 
limit in which we can neglect all non-zero Matsubara frequencies, and (0) reduces to the classical problem in which 
S are commuting vectors of length S . The equations (1q) are analytically solvable and we can obtain a closed form 
expression for the critical temperature at which spin-glass order vanishes: 

^4 - -ij^JS\ (40) 

A second, more sophisticated limit, valid at lower temperatures (well within the spin-glass phase) is when we 
examine uj and T of order JS. It is therefore useful to define the variables uj — uj/{JS) and T = T/{JS), which 
remain of order unity at large S. With this scaling, the integral equations ( pl|) reduce to independent quartic equations 
for each frequency. More precisely, if we make the following Ansatz for the Green's function 

G(a;,r) = jg9i[^,T) + -^g-^i^J) + . . . , (41) 



then to leading order in 1/5*, (31) reduce to 

/oo 
pi(cJ)(iZj 
-OO 



gi{u)-^ =uj---m- 2.gi(cJ) - gi{-u) (42) 

where pi = — Imgi/Tr as usual. Eliminating the frequency —ZJ, we find a quartic equation for gi{uj). We do not 
explicitly display the far more complicated equation for the subleading term g2. 

The solution of the quartic equation for — Q^ is presented on Figure ra together with a numerical solution of the 
full integral equation for S = 5. From the solution of the quartic equation we find that pi vanishes linearly frequency 
at low frequencies; indeed, we find the analytic expansion 

,-, 1 (l + ^) _ , (2-3^)V3 _2 
9i{^) = — ^ ^ — ^+ -^ ^ (43) 

at low frequencies. We expect that the full Green's function in (|4l|) also has a similar low frequency expansion, 
although this has not been proved. It is not difficult to show that the linear low frequency spectral density holds at 
higher orders in the 1/S expansion; moreover, our numerical results, shown in Fig 7|, also clearly indicate a linear 
behavior at small ld. At dominant order in the present large S theory, the spin susceptibility is given by : 

■^^"{u)) = -n dxpi{x)pi{x-uj)+g7:{pi{u;)- pi{-uj))+'Kg'^l3ujS{uj) (44) 







and this is also shown in Figure m. 

It is important to realize that the deceptively simple structure in (|4^) relies on the special value 8 = Qr = 1/a/3 
determined by the entirely different replicon argument in Section III B l[ For arbitrary values of O we either find no 



physically sensible solution of the large S quartic equation (this is the case at the replica symmetric value = 1 
where the spectral density does not satisfy the required positivity criteria) or a solution with a spectral gap. In the 
latter case, the solution for gi is real for small real ZZJ, and there is an onset in the imaginary part ~ (a7 — ZUc)"'^^^ above 
some critical frequency uJc- The solution for = Q^q is of the second type: it has a finite energy gap, but does not 
violate any spectral positivity criteria. 

This subsection has so far identified two distinct large S regimes. In the regime T ~ JS^ we have purely classical 
behavior (the non-zero Matsubara frequencies can be neglected for static properties) and a phase transition at a 
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FIG. 7: &) piiiS) at T K for S — 5 and 9 = Or : the solid line is the numerical solution for the integral equation 
(pll), the dashed line is the solution of the quartic equation (pj|j for gi. b) x" i'^)/^ at T w from the quartic equation. 



critical temperature in ( [40| ) where the spin-glass order vanishes. At lower temperatures, T ~ JS, we are well within 
the spin-glass phase, and the semiclassical dynamics is described by the solution of a quartic equation defined by (H3). 
As we noted in Fig |^, there is a third "quantum" regime at even lower temperatures, T ^ JVS, and this becomes 
evident in a study of the thermodynamic properties presented in the following section. 



D Thermodynamics 



We now turn to the internal energy U and the specific heat C. By computing the average of the Hamiltonian in 
the N —^ oo limit, we find that U is given by : 



U{T) 



2 



/9 



dT[G''''iT)G^\-T)y 



(45) 



Using ( pq ) and the one-step replica symmetry breaking Ansatz in the spin glass phase, we find : 

r/3 



[/(T) 



4^ GirrCi-rrdr 



j2 r^ ~ 



gf{G{-r)-gfdr-:^Pix-l)g^ 



in the paramagnetic phase 



in the spin glass phase 



(46) 



A numerical computation of the internal energy U{T) and the specific heat C{T) is displayed in Figure M for 
= 0_R- The condition for the phase transition between the two phases is that the breakpoint in the Parisi function 
reaches its limiting value x = 1. In this limit the equations determining the parameters in the spin glass phase, ([3]]), 
transform continuously to those for the paramagnet. As the equations are believed to have a unique solution, this 
implies that there is no discontinuity in the internal energy at the transition, indicating its second-order nature. This 
is confirmed by the numerical solution displayed on Figure g. 

Moreover, in the large-S* limit defined above, we can perform a low temperature expansion of the internal energy. 
Inserting (^l|) and (^) into (|4|), and keeping the dominant terms at large S we find 



U{T) 



i + 3e\ g^ - Jh^^Y. [G{i^n)G{-ivr,) + 2G^{iv^) 



^ (^ + Se") +JS(^ + Se") r^5i(*I7„) - JStY, [5i(*^n)5i(-*^n) + 2g?(*I7„)] + 0(J5")(47) 
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FIG. 8: The speciGc heat C[T) and the internal energy U{T) vs. 
Eqs. @) for 5 = 5 and e = 6^, 



the temperature T , from a numerical solution of 



where V„ = Vn/{JS). Clearly, this result indicates that the leading term in U{T) is a temperature- independent 
constant of order JS^, followed by a term of order JS whose coefficient is a function only of {T/JS). Evaluation 
of the latter function at low T for G = Qr yields a curious accident: the gapless structure of the spectral functions 
suggests that the low T expansion should depend only on even powers of T/JS, but it is not difficult to show using 
(Es) that the coefhcient of the term of order SiT/ JS)"^ vanishes. The first non- vanishing, T-dependent term among 
those shown explicitly in ( |47| ) turns out to be order JS{T / JS)"^. To obtain the true low T behavior we need to 
expand (E^ to one higher-order in 1/5*, and this requires use of the second term, g2, in (41). We do not expect any 
cancellation of the term of order [T / JS)"^ at this point, and so the low T expansion for U looks like 

U{T) = t/(0) -I- aS{T/JSf + b{T/JSf + . . . (48) 

Rather than numerically evaluating the values a and &, we will be satisfied by the full numerical solution of (131 
followed by the evaluation of (46). The results are shown in Fig g and are consistent with (|48|). The structure of the 
expansion in (|48| ) suggests that these results are valid for T < J^/S, where the specific heat depends linearly on the 
temperature. Although the present discussion has been carried out for large S, we expect, and this is supported by 
our numerical results, that the linear T dependence of the specific heat holds even for small 5' as T -^ 0. 

In Appendix O we describe the computation of the specific heat of the quantum rotor and Ising spin glasses 
considered in Ref. |3[ As noted in the introduction, these models are simpler because they do not have quantum 
Berry phases in their effective action. Further, at low orders in their Landau theory, the solution for the spin glass 
phase is replica-symmetric. However, understanding the true T ^ behavior requires inclusion of higher-order, 
"dangerously irrelevant" terms which induce replica symmetry breaking; this is carried out in Appendix |Q, and we 
find that these quantum spin glasses also have a linear specific heat at low T. 



IV. CONCLUSION 



We believe that the results of this paper provide a reasonably complete understanding of the infinite-range quantum 
Heisenberg spin glass. While there have been a large number of previous studies of quantum spin glasses of Ising spins 
and rotors (including models with {p > 2)-spin interactions), none of these models contain quantum Berry phases 
in their effective a ction s, as is the case with the Heisenberg model. They have strong consequences: the spin-liquid 
solution of Section [I A and its spectral density IM) are novel properties of the Heisenberg model. There is an intricate 
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interplay in stability between this spin-liquid state and the state with spin-glass order at low T which we have also 
described. At sufficiently low T, the spin-glass order always appears, and we have also described the thermodynamic 
properties of this state. 

An important issue not resolved in our analysis is the origin of the marginal stability criterion in the fluctuation 
eigenvalues in replica space. We imposed this criterion in a rather ad hoc manner, and found that it was the unique 
case under which the quantum excitation spectrum was gapless. Ultimately, the selection criterion for the spin glass 
state has to be a dynamic one, and this requires an analysis of the approach to equilibrium in real-time dynamics. 
Such an analysis was not carried out here, and is an important direction for future research. 

Another interesting open problem is to extend the study of d3) to cases where Jy has a non-zero average value. 
This will allow for ground states with other types of magnetic order, ferromagnetic and antiferromagnetic, and 
their competition with the spin glass state should be of some experimental interest. Interesting transitions in the 
paramagnetic states from the spin liquid state discussed also appear possible. 

We have already mentioned a recent study pq| of the quenching of the spin liquid state by mobile charge carriers 
into a disordered Fermi liquid. Combining this with models just mentioned, with a non-zero average Jij, should lead 
to results of direct physical interest in the heavy fermion and cuprate series of compounds. 
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APPENDIX A: COMPUTATION OF THE SPECTRAL ASYMMETRY 

This appendix is devoted to the derivation of Eqs dlQ). We will consider hereafter the fcrmionic case (the bosonic 
one is very similar). At zero temperature, the number of particles is given by 



qo = z —G''f{Lo)e'^^'=iP -=l9^1nG^He*"-' -zP / — G^ H^^S^ (c.)e'-'^^ (Al) 

where G^ is the Green function with Feynman prescription on the real axis, and the symmetric principal part is 
defined by 



P = hm / + / (A2) 

Using the relation between G^ and the retarded Green function G^ , we find for the first term : 

f°° duj P n+ argG?(0~) -argG?(-oo) f°° duj c ■ n+ 

^P — d^ InG-f (w)e^'"" = ^^ ' ^^ + iP — d^ lnGf(Lj)e''^° (A3) 

J-oo 27r ^ TT J_^ 27r ^ 

The arguments can be extracted from the low-energy and the high-energy behavior of the Green function, which 
leads to argG?(0~) = — 37r/4 — and argG?(— oo) = — tt respectively. The integral on the right of (A3) can be 
easily evaluated : we close the contour of integration, avoiding the singularity aX ld — Q and use the analyticity of the 
retarded Green function in the upper half plane, in which it has no zeros nor poles. We find finally : 

q, = \---^P r '^^ Gf (c.)a^Ef (^)e-o" (A4) 

2 TT J_^ 27r 



The problem is now reduced to the computation of the integral in (A4) as function of 9, which turn out to be the 
most difficult point. An analogous computation was performed in the overscreened regime of a large- A^ description 
of Kondo effects |2^, ^ , but it turns out to be more complex here. Proceeding along the lines of Ref . |2^, ^, we 
note the existence of Luttinger-Ward functional ^lw ~ J dtG^{t)G^{—t) which has two properties : first we have 
I]^{uj) — S^Lw/5G^{u!); second ^lw is invariant in the transformation G{llj) — > G{u! + e). From this, we could 
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naively think that the integral of (A4) vanishes. However, it is not possible to find a regularization for the integral 
for which we could use the invariance of the Luttinger-Ward functional and a more careful analysis shows that : 



,p r ^ r.. 



27r 



G)(uj)d^^{L0) 



sin 20 



(A5) 



To obtain this result, we introduce the following paranietrisation of the singularity at w = 



C_ 



for w > 
for tj < 



(A6) 



The principle of the computation is very simple : we compute explicitly the integral with a regulator 77 > and then 
perform the limit 77 — > 0. Going to the real axis, we find : 



S^(^) 



, , , p{^\)p{0J2)p{uJz) 

dujiduj2auj3 -— 

a;i>o wi<n cji +LJ2 - ^3 - ti^ - iO^sgntJi 

ti'2>0 OU UJ2<0 

a;3<0 iiJ3>Q 



(A7) 



Using the notations a — a -\- ica, b = b + ieb, t^ih = ±0^ and iprjix) — 6 (|a;| — rj) (a are b real and O is the Heaviside 
function), we obtain (using the definition of the principal part (|A2[)) 



Ma,b)=P 



dz 



-00 {z — a) (z — b) 



(A8) 



{a-bf 



In 



{'q + b){ri-a) 



iv-b){ri + a) 
Using the spectral representation for G^ and (A7), we find 



i7r?/',,(6)sgnef, - in ip ^^ (a) sgne a + 



a ^ b \ri — a 77 + a 



j^dwfc p(wo)p(wi)p(w2)p(w3)0j)('^i +u;2 - W3 - ieiSgnwi,a;o -leosgnwo) 

A1UA2 ^^Q 



(A9) 



with an explicit integration over lo with (AS). In this expression, the integration domains are defined as 

u 

[w3<0j [cJ3>0j [cl'3<0J 

Since (j)r]{—a, —6) = 




A2 = 
I ^ I tJ2 < 

W3 < 0^ 
t>ri{a, 6), a simple change of variable leads to : 




(AlO) 



I 



Xi>0 



(p{u}i)p{u}2)p{oJ3)p{^o) - p{^i)p{^2)p{^3)p{i^o)]4>r,{xi + X2 + X3 -iei,-xo + ieo) 
+ (p{uji)p{uj2)piuJ3)p{uJo) - p{^i)p{^2)p{^3)p{^o)]4'r,i^i +^2 + X3 ~iei,Xo ~ leo) 



(All) 



with piu) = p{—Lu). To tak e th e limit 77 — > 0, we use the n ew variables Xi — rjUi and the behavior of p{x) for a; — > 0, 
parametrised according to ( |A6| ). The first integral in (All) vanishes at dominant order in 77 (this term is proportional 
to C\Ct — CtC\ = 0), but the second integral gives : 



ClC--ClC+ 
1 

„.>0 ^/UqUiU2U3 



^^n=i{ui + U2 + U3 - iei.UQ - ieo) 



(A12) 
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Using X — uq, y ~ ui + U2 + U3 and polar coordinates in ^/ul, we find I = 27r(C^C_ — C'^C+)l2 with 



dx 



--Vydy 



{x — y + ieY 



In 


l+xl-y 


l-xl+y 



+ i7r(i/'i(y) - V'i(a;)) 
1 



After an integration by parts on y and using 



^ dx 



In 



V 



T 



1 



1 



I - y + iei 1 + y — iei 



we find 



iP 



d(jj 
2^ 



Gj{uj)d^j:JiLu) 



{clc^ - ClC+) 



Finally, a analogous computation can be performed in the bosonic case, leading in both cases to : 

sin 26* 



iP 



Ztt 



(A13) 



(A14) 



(A15) 



(A16) 



(in this expression, — tt < 9 < tt). These expressions have been shown to agree perfectly with numerical computations 
in imaginary time for the fermionic case and on the real axis at zero temperature in the bosonic case. 

Let us note finally that we can guess the result if we admit a priori that the integral is given by an homogeneous 
polynomial of degree 4 : due to the particle-hole symmetry (in the fermionic case : / <-^ /^ , the result can be expressed 
as a function of C^ — Cl and C'^C- — C'^C+. The first term is rejected since it leads to a singularity at 9 = ±7r/4. 
The proportionality coefficient is fixed by imposing 9 — 7r/4 for qq = 0. 



APPENDIX B: THE MARGINALITY CRITERION 



1 Diagonalization of the fluctuation matrix 

First, we diagonalize the fluctuation matrix M in the replica space defined by Eq. (|33[). A priori M is a. n{n — 
l)/2 X n{n — l)/2 matrix. However, we have taken for gab the simple one step replica symmetry breaking Ansatz on 
gab, 'i-B- the n x n matrix splits into n/m x n/m blocks : gat = 3 if [a/^ = \b/m\, otherwise. Thus M splits 
into n/m identical m(m — l)/2 x m{m ~ l)/2 blocks {Mab,cd does not vanish if and only if all indices are in the same 
block [a/?7ij = [6/toJ = [c/nij = ld/m\). Hence the diagonalization is to be performed only on one block (we set 
1 < a,b,c,d < m), which elements are given by (with a, 6, c, d distinct replica indices) : 



Mab,ab = A= ^pj^g^ 
Mab,ac = B = -9/3V4g5 j 2g 



1 - 3/3^252 ( g2 

e 



e 



pjg 



l3 74 „6 



This matrix has already been diagonalized in Ref. pq and its eigenvalues are given by 

ei=A-2B + C 

62 = A + 2(m-2)B + (TO-2)(m-3)C/2 

63 = A + (m - 4)B - (m - 3)C 



(Bl) 

(B2) 
(B3) 



(B4) 
(B5) 
(B6) 



Moreover, the degeneracies of ei, 62 and 63 are n{m — 3)/2, n/m and n{m — l)/m respectively. Using (Bl) and 
we find finally the result of the text ( p4|) . 

The above calculation has entirely ignored perturbations, SG{t) in the diagonal elements of (p5|). Including these 
greatly complicates the analysis, but a simple observation will suffice for our purposes. Our main attention is on the 
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cross-coupling between SG{t) and the Sgab- A simple consequence of the block-diagonal structure of the gab in the 
mean-field solution is that this cross-coupling has the form 



ST. 



E 



a>b. [a/mj — Ib/m] 



dTSG{T)Sgab- 



(B7) 



Now we can expand Sgab in terms of the eigenvectors associated with ( B4 ) , which were computed in Ref . |36|. The key 
observation is that after the sum over a,b in (B7), the cross-terms corresponding to all the eigenvectors associated 
with ei vanish. Consequently these eigenvectors remain eigenvectors even upon including 6G{t), and the eigenvalue 
ei remains unchanged. A similar argument shows that the eigenvalue 63 also remains unchanged, and only the 
eigenvectors associated with 62 are modified non-trivially by the coupling to 6G{t). 



2 The replicon solution is gapless 

Let us assume that there is no gap in the boson spectral density and more precisely that for small lu : 

G(u;) = -I- -f (a + tb)iu°'e(Lu) + (a' + ib')\Lure(-Lo) + oilton (B8) 

Jg 

where Q is the Heaviside function, a, 6, a' , b' are real constants, and a > 0. Then from ( |31b| ) we obtain for ui > : 

Im (l^{uj)] = J^g^lm OlGiuj) + G{-uj)\ +... (B9) 

(BIO) 



I](w) =c+{d + i{2b - b')J^g^) LO 



where c and d are real constants. The other terms are subdominant in the limit oj — > as can be seen using a spectral 



representation. We then expand (31a) to second order and obtain at first order A 



Jg 



■S- and for the imaginary part 







at second order : 

b~{2b~b')(d^ ^b' -{2b' ~b)Q'^ 
which leads to (6 = 1 is excluded since 6 > and b' < 0) 

Thus the value of Q given by the replicon condition is the only one that leads to a gapless bosonic spectral density. 



(Bll) 
(B12) 



APPENDIX C: FREE ENERGY OF QUANTUM ROTOR AND ISING SPIN GLASSES 

Quantum spin glasses of quantum rotors and Ising spins were studied extensively in Ref. O. However, while the 
paramagnetic phase and the vicinity of the quantum-critical point were fairly completely described, the T — > ther- 
modynamics within the spin glass phases were only studied in the replica-symmetric solution. A proper understanding 
of this low T limit requires consideration of replica symmetry breaking, and we will provide that here. We will find, 
as in the more complex Heisenberg spin model considered in the body of the paper, that the specific heat is linear in 
T at low T. 

As we are restricting our attention to mean-field theory, we can neglect the spatial dependence of all degrees of 
freedom. Further, we will also restrict ourselves to the Ising case, and the generalization to the multi-component rotor 
case is immediate. As discussed in Ref. |l3|, the effective action of the quantum Ising spin glass is expressed in terms 
of the order parameter functional 



where cr° is the Ising spin in replica a. The important low-order terms in the free energy density are 

d d 



(CI) 



T 



U-^ 



dri dT2 



Q""(ri,T2) 



Ti=T2=T 



^ J dT,dT2dT-i Y, ^"''(Tl, T2)g^^(r2, T3)Q^"(r3, Tl) 
abc 



I J dTY,Q^^{r,T)Q'^''{T.T)-y- j dT^dT2Y,[Q^\T^,T2)\\ 



(C2) 
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Here r is the parameter which tunes across the spin glass transition, and k, y measure the strength of various non- 
Hnearities. The analysis of thermodynamic properties in the spin-glass phase in Ref. |l3| was carried out with a 
vanishing coefficient of the quartic term, y = 0: in this case the order parameter has replica symmetry, and it was 
found that the specific heat ^ T'^ as T — > 0. Here, we will extend the solution to small y ^ 0, and show that the 
solution with broken replica symmetry has a linear specific heat. 

Time-translational symmetry requires that the mean-field solution take the form 



Q"''(Ti,r2)-i^Q''''(*z.„)e'''" 



(ri-T2) 



(C3) 



As in (|25|), we choose the following Ansatz for Q 



{ivn) 



D{iUn) + PqEA 

(3qab 



a = b 



(C4) 



where the off-diagonal terms, qab, are time-independent and characterized by the Parisi function g(u), and q{\) = qsA- 
We have included an additive factor of f3qEA in the diagonal term for convenience, and without loss of generality: as 
in the discussion below (Eq), we will find that this ensures that at T = the solution for D(t) vanishes as r ^ cx3. 
Also, the diagonal components qaa do not appear in the above, and we are therefore free to choose them as g"° = 0. 
Here, and in the remainder of this appendix we are assuming that r is sufficiently negative so that the system has a 
spin glass ground state; for larger r, the ground state is a paramagnet |13 with qEA = qab = whose properties are 
not addressed here. 

We now need to insert (C4) into (C2) and find the saddle-point with respect to variations in the functions q{u) and 
D{iVn). This is, in principle, a straightforward exercise, but the computations are somewhat lengthy. 

We first identify just the terms that depend upon q°'^\ these have the form 



T=-R,Trq--^Trq^-^Y.<iab 



(C5) 



ab 



where 



Ri = /3k{D{0) 

i?2 = K/3^ 

i?3 = Py- 



(iq 



EA) 



(C6) 



We first address the problem of determining the saddle-point of (C5) with respect to variations in q{u). Fortunately, 

|, and we can directly borrow the 
4i?i_R3/i?2)"^^^, where it saturates 



this problem has been completely solved in the classical spin glass literature |g 
results: the function q{u) increases linearly as a function of u for < m < 1 — (1 
until u = 1 at the constant value 



qEA 



i?2 



im 



4i?ii?3 



,1/2 



2i?3 



(C7) 



Combining (C7) with (C5), we obtain the simple result 



2 
lEA 



-nD{Q)/y 



(C8) 



Next, we consider the variation of T in (|C2| ) with respect to D{ivn)- This is most easily done for Vn 7^ 0, for which 
we obtain the following saddle-point equation 






^I]^K)+9£A 



"^Vq^sADi-iVn) 



'iyq 



EA 



22/ 



^ ^ D{lv:,)D{-lU,, - IV\,) - ^ E ^{^ODiiyDBi-lVr. - ll^'n - ^O = 0. 



(C9) 



Upon consideration of the saddle point equation for -D(O) one initially finds a number of additional term associated 
with the coupling of i^(0) to the qab- However, our parameterization in ( p4|) w as chosen judiciously, and has the 
feature that all these additional terms vanish upon using (pa); so, the result (|C^) applies also for v^ = 0. 
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Let us also note the complete expression for the free energy density, obtained by inserting ( |C4| ) and (C9) into (|C 

-\ 2 



T/n = 



K 5k (3k 



i-J:i-n + r)Di^.,,)-^^J2D^i^., 



3(3 



■o^Di'i-'^n)+qEA 



P 



yiEA 

p 



y^D{ii'n)D{-ii'n) 



2yqEA 

3/32 



J2 D{iv^)D{iv'^)D{~ 



V 
6/33 



^ D{iVr.)D{iv:,)D{iv';^)D{-iUr, - iv'^ - ivl). 



(CIO) 



We are now left with the task of solving the saddle-point equations ( |C8|) and (C£) for qEA and D[ivr()^ and inserting 
the result in (CIC). This is clearly a daunting task, and we will be satisfied in de scribing the T — > limit to first 
order in y. This is similar in spirit to the large S expansion of Sections [II C S , IIID| , and we expect that higher-order 
corrections in y will not modify the nature of the low T limit. 

First, we consider the case y — 0- Here a complete analytical solution is possible, and was presented in Ref. |l^. We 
have at y = 0: 

qea = i5-yii^»i — 

(3k ^ — ' KU 



T''{T)/n 



3(3k'' 



■El 



2k^u 



^T^{0)/n 



4^ 
AbK^ 



(Cll) 



We observe that the free energy density behaves as T"', while the specific heat ~ T'^ as T ^ 0. Notice also that 
positivity of q%j^ requires an upper bound on r, which we have assumed to hold. 

Before considering explicit corrections in powers of y, we make an observation that is valid to all orders in y. The 
solution for D(ivn) in (Cll), when analytically continued to real frequencies, w, has an imaginary part which vanishes 
linearly in uj at small uj. We now show that this conclusion holds to all orders in y; the constraint ( |Cq ) will play a 
key role in establishing this result. Let us write D{uj) = D{0) + iDiUj + . . . for small lo, where D{Q) and Di are some 



real constants. Inserting this in (C£) and evaluating it at T = for small w, we note that the last two terms in (C9) 
have imaginary parts which vanish as oj'^ and lo^ . Keeping only the leading to dependence of the imaginary part, we 
obtain the simple expression 



-2iKD{{))DiUj - 2yq%^iDiLU = 0. 



(C12) 



From (C8) we see that this condition is satisfied, and so Di can be non-zero. 

Now, we consider explicit first-order corrections in y: we w ill sec tha t this leads to terms in the thermodynamics 
which vanish more slowly as T ^ 0. We can easily use (C8) and ( |C9| ) to determine the corrections to D{ivn) and 
qEA to linear order in y; however, these are not needed here as the shift in the free energy due to such corrections 
will only appear at order y^, because t he fr ee ener gy is at a saddle point. Indeed, to obtain the free energy correct 
to first order in y, we need only insert (Cll) into ( |C10 ). It is then quite easy to see that the free energy will have a 
term of order yT"^ , and that the coefficient of this term will be non- universal and depend ent u pon the nature of the 
high energy cutoff. The required term comes from the T dependence of q%j^, which from (Cll) is seen to be 



4a(T) = 4a(0) 



3k 



(C13) 



Upon inserting (C13) and ( pll ) into (CIC) we will now obtain numerous terms in which the above T^ term multiplies 
T-independent, cutoff-dependent terms coming from the upper bounds in the summations over the D{ivn). As the 
relative values of these contributions will depend upon the nature of the cutoff, there is no general reason for them to 
cancel against each other. Hence we obtain a T^ contribution to T and a linear T term in the low T specific heat. 
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